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Current and noise expressions for radio-frequency single-electron transistors 
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We derive self-consistent expressions of current and noise for single-electron transistors driven by 
time-dependent perturbations. We take into account effects of the electrical environment, higher- 
order co-tunneling, and time-dependent perturbations under the two-charged state approximation 
using the Schwinger-Kedysh approach combined with the generating functional technique. For 
a given generating functional, we derive exact expressions for tunneling currents and noises and 
present the forms in terms of transport coefficients. It is also shown that in the adiabatic limit 
our results encompass previous formulas. In order to reveal effects missing in static cases, we 
apply the derived results to simulate realized radio-frequency single-electron transistor. It is found 
that photon-assisted tunneling affects largely the performance of the single-electron transistor by 
enhancing both responses to gate charges and current noises. On various tunneling resistances and 
frequencies of microwaves, the dependence of the charge sensitivity is also discussed. 
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I. INTRODUCTION 

For a metallic island small enough to give large charg- 
ing energy exceeding the temperature, quantum trans- 
port through it shows remarkable features due to strongly 
correlated electrons. A single-electron transistor (SET) 
is popular geometry of studying transport through the 
metallic island, in which the island is coupled to two 
large reservoirs (source and drain) via tunnelingiunctions 
and to another reservoir capacitively (gate). 0Xl El ELS 
The most important principle of operating SET is the 
Coulomb-blockade effects in the island. In an usual sit- 
uation particles cannot tunnel into the island due to the 
Coulomb energy. However, when the Coulomb energy is 
reduced by a gate bias, transport through the island is 
possible even with a small bias between the source and 
drain. Consequently, tunneling currents show a series of 
peaks, called the Coulomb-blockade peaks, as a function 
of gate voltage. 

From fundamental and applied point of views, the 
shape of Coulomb-blockade peaks has been attracted 
much attention and widely studied in static bias 
conditions. p The appearance of the Coulomb- blockade 
peaks depends on two conditions. The first is much larger 
charging energy than the temperature to blockade ther- 
mal excitation. The second is thick tunneling barriers 
guaranteeing a large dwelling time for electrons to re- 
solve the charging energy in the island. Then, from the 
energy uncertainty principle, the latter is approximately 
fulfilled under the condition that parallel resistance Rt 
of the barriers is much larger than the resistance quan- 
tum R K = h/e 2 , i.e., a = Rk/{^ 2 Rt) < 1. It is 
now well known that for large tunneling resistance Rt, 
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transport in a SET is achieved by the sequences of un- 
corrected tunneling processes and the smearing of peaks 
is dominated by the temperature. However, for relatively 
small tunneling resistance (ao < 1) or very low tempera- 
ture, additional quantum processes such as higher-order 
co-tunneling contribute the peaks and renormalizes even 
its positions as well as associated quantum fluctuation 
is responsible for the broadening of the peaks. 0, Q] It 
is also found that these higher-order co-tunneling effects 
are manifested to noises of the system. 0, H, On the 
other hand, in concerning about applications of SETs to 
an electrometer, a slope of a Coulomb-blockade peak with 
respect to gate voltages is an important factor. Since the 
distance between peaks is equal to the change of an el- 
ementary charge on the gate, a large slope of the peaks 
means high sensitivity to a fraction of the charges. So, 
it is widely believed that a SET is a prime candidate 
for reading out the final state of a qubit in a solid-state 
quantum computer. 0, 0, E, 0] 

New theoretical interests in transport properties 
through the strongly correlated systems are emerging 
together with the experimental success in driving them 
by microwaves (radio-frequency waves) , which are called 
radio- frequency single-electron transistors (rf SETs1.[l5j 
In such a system, microwaves are delivered via a LC- 
resonant circuit to excite particles to overcome the 
Coulomb energy. As a consequence, SETs can oper- 
ate in a high-frequency domain and practically pro- 
vides advantage of a large bandwidth as an electrome- 
ter, which allows to measure the rapid variation of gate 
charges. [H E E [H 03 Theoretically, the interplay 
of electronic transport and excitations by microwaves is 
a particular interest because high-frequency perturba- 
tions are expected to yield a new non-equilibrium situ- 
ation resulted from additional phase variation in energy 
states. Such a time-dependent situation is usually di- 
vided into classical and quantum regimes. In the classical 
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regime (or adiabatic regime) energies excited by time- 
dependent perturbation appear to be continuous while 
in the quantum regime (we will also refer to this as non- 
adiabatic regime) discrete photon energies become ob- 
servable and particles can emit or absorb photons when 
they tunnel from an initial state on one side of the barrier 
to a final state on th e oppos ite side, called as photon- 
assisted tunneling. [20 II21L Ej I23 ] 

So, in order to understand transport properties of rf 
SETs, one may need generic theoretical considerations 
including higher-order co-tunneling processes as well as 
sequential tunneling, even in the quantum regime of time- 
dependent perturbations. Actually, according to the re- 
cent experiments of rf SETs. [lallallall7lll8| tunneling 
resistances range from ao = 10~ 5 to 3 x 10~ 2 , implying 
feasible co-tunneling processes for large values of ao- Fre- 
quencies of microwaves were used from 0.3 to 1.7 GHz, 
which correspond to several /ieV of photon energies, com- 
parable or larger than thermal energy in the experiments. 
For these values of frequencies, it has been not known 
for the system to be driven in the classical or quantum 
regimes of time-dependent perturbations. So, in general 
it is necessary to solve problems in the quantum regime 
for more rigorous understanding of transport in rf SETs. 

Additionally, tunneling in a rf SET may be dissipa- 
tive due to a LC resonant circuit. Since a microwave is 
delivered via a coaxial cable with 5017 impedance much 
smaller than the resistance quantum, effects of the elec- 
trical environment is usually ignored. This is the case for 
a LC resonant circuit with a low quality factor, however, 
with a high-quality factor, energy states of the electri- 
cal environment become long-lived because it becomes 
similar to a simple-harmonic oscillator. 0, IE EH HE US 
Then, during tunneling, particles may emit or absorb en- 
ergy quanta equal to resonant energy of the harmonic 
oscillator to the environment. 

In this work, we develop the formalism that is capa- 
ble of treating all above theoretical considerations; the 
effects of higher-order co-tunneling, non-adiabatic time- 
dependent perturbations, and the electrical environments 
on operations of rf SETs. Our work is a generaliza- 
tion of several earlier works which address the effects 
partially, neglecting time-dependent perturbation. [Toj 
higher-order co-tunneling, [^J and electron-electron 
interaction. [23,|2^| However, since we use a two-charged- 
state model in a metallic island assuming large charging 
energy, along this aspect, our work is more restricted 
than Ref . [3, 0, |2^| . In solving the problem, we use 
the Schwinger-Keldysh approach combined with a gen- 
erating functional [l0|, |2jJ |30j, |3l| where pseudo-spins of 
two-charged states are treated with the drone-fermion 
mapping. Since this approach includes any higher or- 
der moment of diagrams systematically, it is one of the 
well-suited methods for transport through strongly cor- 
related system as indicated in Ref. [lj}. From a gener- 
ating functional summed diagrammatically, observables 
of the system are obtained by functional derivatives with 
respect to external perturbations. This is another ad- 



vantageous point of this approach because higher-order 
moments such as noise are easily calculated and expres- 
sions for observables are consistent with each other in a 
sense that they are derived from the same order of dia- 
grams. We express the electrical environment in terms 
of infinite number of driven harmonic oscillator follow- 
ing Caldeira and Legget 32] where external alternating 
voltages are treated with classical fields. Based on the 
unitary transformation which leaves the electrical envi- 
ronment in a stationary situation, we incorporate equilib- 
rium fluctuation of the environment into the generating 
functional and derive environment- and time-dependent 
self-energies by counting dominated diagrams. Results 
for currents and noises are expressed in terms of trans- 
port coefficients. In cases of time-dependent perturba- 
tions, due to the displacement component currents are 
found to depend on an additional transport coefficient, 
leading to a generalization of the Landauer formula[33j 
and noises also has its contributions. 

The paper is organized as follows. We first describe the 
model of calculations in Section II. By expressing the dis- 
sipative environment in terms of driven harmonic oscil- 
lators, we give the Hamiltonian depending on tunneling 
currents. In Section III, we calculate an approximated 
generating functional based on the Schwinger-Keldysh 
approach and discuss several approximations in deriv- 
ing it. For a given generating functional we show ex- 
act expressions for currents and noises in Section IV and 
rewrite them in terms of transport coefficients. In Section 
V, we use our formalism to simulate a rf SET numerically 
and emphasize different points from static calculations. 
Finally, we summarize our results in Section VI. 



II. MODEL OF CALCULATIONS 

A. Hamiltonian 

To formulate the problem of a rf SET, we begin 
with general circuital geometry of Fig. ^ where time- 
dependent external sources, Vo(t) and Vc(t), are sup- 
plied via dissipative elements of impedances Zd(uS) and 
Zq{lo). In this section we do not specify detailed forms 
for Zd(lo) and Zq{uj) bearing in mind the application 
of our formalism to othe r sy stems concerning effects of 
dissipative environments. |34 l35l l36| As a typical model 
of a SET, a small island is coupled via tunneling barriers 
to two leads, source and drain, and also capacitively to 
source, drain, and gates with capacitances of Cg, Co, 
and Cg, respectively. We assume that the small island 
is a metallic one, i.e., there are many energy levels with 
negligibly small level spacing and also many particles oc- 
cupied to them. In such a metallic island, one can treat a 
excess charge of Q confined in it as a independent variable 
from those of quasiparticles in a good approximation, and 
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FIG. 1: A typical drawing of the single-electron transistor is 
shown where time-dependent voltages are applied to a quan- 
tum dot via possible dissipative elements connected to drain 
and gate electrodes, respectively. 
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where a n ki (Q> n ki) are the annihilation (creation) oper- 
ators for quasiparticles with energy eki in the island 
and the index n describes the transverse channels in- 
cluding spin. The second term is a Coulomb-blockade 
model of the electron-electron interaction with Cs = 
Cd + Cs + Cq- Further simplification of the Coulomb in- 
teraction term can be made if one uses a two-state model 
for excess charges. Assuming the small island enough for 
charging energy Q 2 /2Cs to be the largest energy scale 
in the problem, it is sufficient to consider two number 
of charged states, say, | 0) and | 1). Then, the charge 
operator Q becomes Q = e | | and satisfies Q 2 = Q 
(throughout the work e is the proton charge) . If we adopt 
spinor notation, then the Hamiltonian is further written 
as, 



Hi - ^ e k ia nkI a n ki 



1 



(2) 



nk 



where a z is the effective spin- 1/2 operator and A = 
Ec(l — 2qo/ e ) is the energy difference between the two 
charge states together with the charging energy E G — 
e 2 /2Cs. Here, we anticipate An, which depends on a 
static component of a gate voltage Vq through a charge 
of q = C G VS. 

As for the remaining parts of the system, the Hamilto- 
nian can be found by separating it into terms depending 
on macroscopic and microscopic variables in a similar 
manner to Ref. j2^. Then, the total Hamiltonian may 
be written as the sum of the unperturbed part and the 
tunneling part H.' T ; 
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(V D ,V G ) + H' T 
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Here, TLi, Hi ea d, and Hrlc(Vd, V g ) represent the un- 
perturbed part of our system, which are shown in Fig. 

121 The Hamiltonian Hi ead = J2nk,£=s,D e kta n u a nU de- 
scribes non-interacting electrons in the source and drain 
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FIG. 2: Microscopic and macroscopic parts of the single- 
electron transistor in Fig. are shown. Here, we define 
Z' g (lu) = Z g (lo)(C d + C S ) 2 /C 2 D and V G (t) = V G (t)(C D + 
C S )/C D . 



with their annihilation (a n ke) and creation (a nkg ) oper- 
ators while the Hamiltonian Hrlc(Vd,Vc) governs the 
electrical environment which corresponds to a lumped 
circuit of Fig. [3 Actually, the lumped circuit is de- 
signed to exhibit the same dynamical behavior as that 
of Fig. 1 if it were not for tunneling and thus tunneling 
barriers work as just capacitances. For this, dynamical 
variables in the two circuit are related to each other as, 
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where Qi (I = D,S, G) are excess charges on the ca- 
pacitors in Fig. 1 and 4>i corresponding phases which 
are related to the potential difference vi(t) via the re- 
lation of 4>t(i) = evi(t)/h. Here, the capacitance C\ 
and Ci are defined by Ci = CdCs/(Cd + Cs) and 
C 2 = C 2 n C G /Cv(C D + C s ). Following Caldeira and 
Leggett|32|, one can express the circuit of Fig. 2 in terms 
of many coupled harmonic oscillators each of which is 
quantized under the commutation relation of [4>j, Qj] — 
ie. So, considering classical fields of voltage sources ap- 
plied to the lumped circuit, the electrical environment is 
described by a set of driven and coupled harmonic oscil- 
lators. 

The tunneling part of the Hamiltonian may be given 
as,0 



n'At) = E E [ : 
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where T^i, denotes an element of tunneling matrix be- 
tween a state | nk) in the lead £ and a single particle 
state | nk') in the island, and usually approximated as 
Tui, ~ T l independently of energy levels. Here, the oper- 
ators of cr+ and e 1 ^ 1 are inserted for the increase of excess 
charges in the island and the lead £, respectively. From 
the commutation relation the operator e^ e can be shown 
to increase excess charges by the elementary charge e in 
the lead I for every tunneling event. 



B. Current and self-consistent calculations 

The current Ie(t) flowing in each lead of Fig. ^is de- 
fined to be positive if it flows into the island, so that 
some of charges carrying the current are used to increase 
charges on the capacitor connected to the lead while the 
others tunnel into the metallic island. The former repre- 
sents the displacement current if it) which is equal to a 
time-derivative of the averaged charge as, 



(J) 



where Qe.it) is the Heisenberg representation of the 
charge operator Qi and (• • • )o stands for the ensemble 
average. Whereas, the latter is the tunneling current 
I\it) which is equal to a time-derivative of the averaged 
particle number as, 



e^W))o 
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where Afe(t) is the Heisenberg representation of the num- 
ber operator Mi — J2 n k a nke anke a * the lead I (we assume 
particles as electrons). In other words, the current Iiit) 
can be regarded as the sum of the displacement current 
and the particle current which are contributed from the 
macroscopic and microscopic system, respectively, 



hit) = Ifit) + Ij(t) i£ = D,S,G), 



(9) 



with Iq = 0. These currents automatically satisfy the 
current conservation relation of Id it) + Is it) + Iq it) = 
despite of time-dependent perturbations, as emphasized 
in Ref. [4l|. This is another consequence of the charge 
conservation that a Gaussian surface enclosing three ca- 
pacitors defining the island always contains zero total 
charges as inferred from Eq. QJ. Then, from the conti- 
nuity equation for charges, the sum of the currents should 
be zero, alternatively, d(Qit))dt = I D it) + Ig it) implying 
the fact that the increase of charges in the island is en- 
abled only by tunneling processes. From the Heisenberg 
equation of motion for Qit), one can show that this is 
the case for our system. 

As for the displacement currents, it is possible to ob- 
tain further analytic forms because the macroscopic sys- 
tem of Hrlc consists of linear elements. By viewing the 
tunneling currents as another external sources as well as 



voltages of Vb(£) and Vc(i), from the Heisenberg equa- 
tions of motion for Qi(t) it is straightforward to show 
that, 
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where each current is expressed in its Fourier component 
defined by lit) = J e lu)t Iiu))du). Here, the voltages V D 



and Vn are given by 
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describing effective voltage lowering by tunneling cur- 
rents. 

With the above expressions for the displacement cur- 
rents, the problem is now reduced to obtaining the tun- 
neling currents of Eq. (jSJ). For this, it is convenient 
to transform the system in such a way that the electri- 
cal environment of Hrlc leaves in a stationary condi- 
tion. Then, under such a situation, harmonic oscillators 
in Hrlc are expected to vibrate about their station- 
ary positions. This in turn is helpful to assume them 
in equilibrium independently of tunneling events even 
though noises from tunneling may modify their fluctu- 
ation slightly. According to Ref. |23j|. it is possible to 
find the unitary transformation which rotates the system 
by voltages of SVd and 8Vq- By these voltages we mean 
that the system is rotated to have the lowered voltages in 
the macroscopic system by the amounts while correspond 
phases in the tunneling Hamiltonian appear additionally. 
Namely, the rotated Hamiltonian becomes, 

Hr = Ho + Ht, 

Ho = Hi ead + Hi + H RLC iV D - SV D , V G - SV a ), 

l=D,Snk,k' 



Here, additional terms of piit)i£ = D, S) in the tunneling 
Hamiltonian actually describe the external phase differ- 
ence forced by the voltages SVd and SVg in the absence of 
tunneling. In other words, it is related to the correspond- 
ing potential difference v^it) across the tunneling barrier 

from the island to the lead I via pe.it) = e/h d,Tv\ir). 
The potential difference v^it) is given by, 
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where the impedance matrix is defined as, 



Z(«) = 
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with Z' G = Zq{\ + Cs/Cd) 2 - Since the number opera- 
tor Mi is found to be invariant under the rotation, the 
transformation by (SVn, SVg) = (Vd —V^,Vg — Vq) may 
give rise to the simplest situation in calculating the tun- 
neling currents. In this case, the electrical environment 
is in stationary conditions as implied from Eq. (|1 1|) . so 
that the first moment of its dynamics does not influence 
the tunneling currents, but the second moment, at least, 
starts to work. 

Instead of this benefit, the rotated Hamiltonian now 
depends on the tunneling currents. This implies that 
observables from the Hamiltonian also depends on the 
tunneling currents unless both Vq and V G in Eq. I|13|) are 
zero. Especially, in cases of the tunneling currents this 
requires self-consistent calculations to obtain it. Detailed 
forms in the rotated frame become, from the Heisenberg 
equation of motion for Mi, 

I\{t) = (l e (t)) , 

I e {t) = W{t, -oo)J t {t)U{t, -co), (15) 

showing self-consistent behavior. Here, the time- 
evolution operator U (t, to) is defined in the rotated frame 



where (3 = l/ksT is inverse thermal energy and e n i = 
e n — pP„ with uP e an equilibrium chemical potential at a 
lead ijni 0E3 

As for the bosonic system of Ti-RLC^yh^ ^g)> ^ i s now 
under a stationary condition; (<pi(t))o — (Qi{t))o = 
0. Then, its dynamical behavior is characterized 
by time-correlation functions between variables such 
as ((f>i(0)(f>ii(t))o. In thermal equilibrium, the time- 
correlation functions are easily evaluated by exploiting 
the fluctuation-dissipation theorem, and results are, (3, 

El 



du>?RZ u ,(uj) e l < 
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(19) 



where Zw (lu) is related to each component of the 
impedance matrix Z(oj) in such a way of Zdd = Zn, 
Zds = Z12, etc. 



III. STATISTICAL AVERAGES OF 
OPERATORS 

A. Generating functional 



U{t, t ) = T exp (i jf dTW fl (r)) , 



(16) 



where T is the time-order operator. Once the tunnel- 
ing currents are calculated self-consistently, other observ- 
ables depend on them explicitly. For example, current 
noises are calculated as, 



Sit (t, t ) = ({Sli(t),Sli,(t )})o 



(17) 



which is defined by the auto-correlation function of the 
current fluctuation operator, 5Ii(t) = Ti(t) — (Xi(t))o- 



C. Equilibrium properties of reservoirs 

As shown in the unperturbed Hamiltonian, our sys- 
tem consists of three fermionic (Tie, & = D, S, I) and one 
bosonic (TLrlc) systems. We assume these systems as in 
equilibrium independently of tunneling because, due to 
large degrees of freedom, effects of tunneling on their fluc- 
tuation are expected to be negligible. Then, for the non- 
interacting fermionic systems, their dynamics are char- 
acterized by single-particle Green's functions. For a | n) 
state with its energy e„ their explicit forms are as follows; 

S^(M') = ^([M*)»aL(*')])o 

ih 2 
g*(t,t') = l0(i-i')<K(t),at(t')})o 



l-eu~t')e e ^ t - t '^ lh 

ih 



To evaluate the ensemble averages for the tunneling 
currents and noise of Eqs. fTK|l and (fTTj) . we use the 
Schwinger-Keldysh approach combined with the gener- 
ating functional technique. [ToL l30l l3l| We are interested 
in calculating the expectation value of O defined by, 

(O(t)) = Tr{p U\t, -co) O U(t, -co)} (20) 

where po is the grand canonical density operator describ- 
ing the system in equilibrium at t 



-co as, 



po = exp{-/3(Wc 



e=D,s,i 



(21) 



and Z is the equilibrium partition function. In order to 
evaluate the expectation value of Eq. I|2(J|I , we introduce a 
generating functional W — —ih \nZ as an extension of the 
Gibbs free energy. Here, Z is the generalized partition 
function defined as, 



Tr{p U_ ( — oo, oo)U + (oo, -co)} 



(22) 



«tf(t,f) = g**(t',t) 



(18) 



where, by different subscripts in the forward W+(co, -co) 
and backward JA\_ ( — co , oo) evolution operators, we mean 
different external fields applied along each time-branch, 
respectively. Such different fields are usually coupled to 
the conjugate variable of O in the Hamiltonian and, at 
the final stage of calculations, are set to be identical. A 
more compact form of the partition function is obtained 
if we view the inverse temperature as an imaginary time 
like, 

r 1 r—°°-iftP 
Po = expj— / dT(H - ^2 PiMi)j 

1 e=D,s,i 
= U T (-oo -ihP,-oc). (23) 
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Here, in accordance with Eq. (|16|l . we define this den- 
sity operator with the evolution operator U T . Then, the 
partition function becomes, 

Z = Tr{U T (~oo — ih(3,—oo)U^(—oo,oo)U + (oo,—oo)} 
= Tt{U c }, (24) 

and can be interpreted as describing successive evolutions 
of states along C+-, C_-, and C r -time branches as shown 
in Fig. 3. We designate these time branches simply a 
closed-time path C bearing in mind that, along each time 
branch, different Hamiltonians govern the evolution of 
states; that is, Ho + Ht along the C±-time branches 
with different fields and Ho — J2e VtNi along the C r -time 
branch. 

Once the generating functional W is given, the ensem- 
ble average of Eq. (|20|) is obtained through functional 
derivatives of W with respect to external fields. In cases 
of the tunneling currents and their noises, the phases 
Pt.it) are assumed to be different along each time branch 
such as, 



pi(t) = Pt(t)- 

Pe(t) = Pi(t)- 



2 

A Pe (t) 



for forward 

for backward (25) 



where Api(t) is fictitious and will be zero at the final 
stage. Then, the functional derivatives of the evolution 
operators with respect to the fictitious field give, 



<5W+(oo, — oo) 



SA Pe (t) 
SU_(oo, — oo) 



2e 



U(oo,t)J e (t)U(t, -oo) 



Ap e =0 



= —U(-oo,t)J e (t)U(t,oo)(26) 
2e 



and, using these relations, it is straightforward to show 
the tunneling current of Eq. (|15J) to be, 



sw 



h 5Ape{t) 



(27) 
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In a similar way, the current noises are obtained by the 
second derivatives, 



Sw (t,to) : 



2e* r 

~ih\ 



6 2 W 



S 2 W 
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ih I Sp e (t)Spj,(t ) Spf(t)5p e , {t ) 
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6Apt(t)SAp e (t ) 4 Spt(t)S Pl > (t 



(28) 



accompanied with Api(t) — Api(to) = finally. Due to 
the normalization of the partition function Z(Api = 0) = 
1, the second term in the last line of the above equation is 
equal to zero. Nevertheless, we keep this term to circum- 
vent the uncertainty related to the order of operators. 

In order to get the averaged charge in the island and 
its fluctuation, we add a fictitious field to the excitation 



FIG. 3: Time contour to evaluate the partition function 



energy in such a way of Ao ± Ah(t)/2. Then, using sim- 
ilar procedure for the tunneling currents, one gets the 
ensemble average of a charge operator Q as, 



(Q(t))o = - 



zSW 



SAh(t) 

while its fluctuation is given by, 
S Q (t,t ) = ({SQ(t),5Q(t )}) 



(29) 
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Ah=0 



Here, we define the charge fluctuation operators as 
5Q(t) = Q(t) - (Q(t)) . 



B. Evaluation of Generating functional 

We evaluate the generating functional W with the 
coherent-state functional integral method defined on the 
closed time-path of Fig. 3. As an usual path integral, 
the evolution operator Uc{t) is divided into a number of 
infinitesimal steps on the time-path C, and then a resolu- 
tion of the identity is inserted at every time step. In the 
coherent state functional integral differently from usual 
ones, the identity operator is given in terms of eigenfunc- 
tions of annihilation operators instead of coordinate- and 
phase-state basis functions. As a consequence, the eval- 
uation of the partition function Z is reduced to the path 
integral in coherent-state variables over the exponential 
of the action along the time contour C.|2^| In our case, 
the result is summarized as 



Z = Z (((, 



.-As,, 



lUi, 



(31) 



in other words, the exponential of the action Si„ t is av- 
eraged over all reservoirs. The action Si nt describes the 
interaction among reservoirs, which is given by Si„t — 
§ c cLtHt{t) on the closed time-path C. Here, Hr(t) 
is a counterpart of the tunneling Hamiltonian for each 
time branches, and is obtained by replacing all operators 
with their coherent-state variables (complex or Grass- 
man numbers) while Hx{t) = in the time branch C T . 
As for spins, in order to utilize the coherent-state repre- 
sentation, we map the effective spin-I/2 operators onto 
two fermion operators c and d (dronc-fermion represen- 
tation), i.e., (T + = c^(d) + d) and a z = 2c t c — 1.|10| 

By the bracket notation, we mean an average weighted 
with the action of a certain reservoir. The detailed form 
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is defined, for instance over Hi, 



as 



(O) 



7~tlead 



■7 lea d 



V[a r u{t),a* rU {t)]e-^ s ^ O (32) 



where Z ead is a normalization factor implying (1)-H icad = 
1 and equal to the equilibrium partition function of TLi ea d- 
Here, {a n ki(t) , a^ ke (t)} are Grassmann variables asso- 
ciated with their fermionic operators and they satisfy 
anti-periodic boundary conditions; a n ke(—co £ C+) = 
— anki(— oo — ih{3 £ C T ). The unperturbed action of the 
leads Si ea d is given as, 



S, 



lead 



h I dT^^2a* nke (T)^-a kn i(T)-Hi ead {T)Y33) 



using the trajectory notation, in which the function 
Hiead{t) represents Hi ead and Hi ea d ~ X^aWV* in tlic 
time branches C± and C T respectively. 

Since the unperturbed actions are quadratic, the ther- 
mal average over all reservoirs of Eq. ffiljl is reduced 
to Gaussian times polynomial integrals if one expands 
e -S int /ih m £ Q a p 0wer series. Then, using a standard 
procedure of a Gaussian integral, each term in the series 
can be evaluated analytically. Firstly the result over the 
reservoirs (£ — D , S ', I) is summarized by the appearance 
of a particle-hole Green's function be(t,t') in the action 
Si n t. The Green's function b^(t,t') has a form of, 

bi(t,t ) = -ihN ch J2\T e \ 2 9ki{t,t )g H (to,t) 

k 

e ipt(€)-ip t (t a ) ) £ = D,S (34) 

Here, the free-particle Green's functions <?feg(i, to) repre- 
sents the inverse function of their free actions and their 
physical representations are equal to Eq. i|18[) which can 
be obtained by the Keldysh rotation. Then, the function 
bg(t,t') represents just a single particle or hole creation 
in the island by tunneling through a barrier I. Actually, 
in obtaining bt(t, t') we take into account only sequential 
processes of single particle or hole creation; one can see 
it if expanding the exponential of Si n t into a series. This 
is a good approximation for a large number of transverse 
channels N c h, so called, the wide junction limit because 
the sequential particle or hole creation is dominated to 
the simultaneous creation of both, at least, by iV c fo.0 

These sequential processes are correlated by further 
evaluation of the partition function over the c— and d— 
fields. The result reads, 



Z = Z ( exp [ - £ -Tr{(g c f )"}! ) 

exp{-ih / dtdt r](t)g d (t,to)r](t )} 

Jc n= o 



(35) 



where the simplified notation of Tr{g c -^-B-^-} stands for 



S S f $ s 

Tr{g c — B — }= / dMt g c (t,to)^ TrT B(t ,t) TT ^(m) 



and B is the particle-hole Green's function combined 
with effects of the electrical environment by 



BO, to) 



E 



e- lMt) b e (t,to)e lMto \ 



(37) 



Here, g c and g d are equilibrium Green's functions for the 
c— and d— fields with their eigenenergies Ao and zero, 
respectively. The Grassmann field rj(t) is introduced as a 
linear source to the d-field, which gives an additional term 
of — ih J c dt[rj(t)d(t) — d* (t)r](t)] into the unperturbed ac- 
tion. 

Finally, we evaluate the partition function over the 
electrical environment. To do this, we transform bi- 
linearly coupled harmonic oscillators in H-rlc into inde- 
pendent ones by an unitary transformation. Then, since 
each simple harmonic oscillator is assumed to be in equi- 
librium, one can exploit the Wick's theorem to show that 
its thermal average can be expressed in terms of, at most, 
the second order phase correlation. To each term of the 
series in Eq. (|35|) . the application of the theorem leads 
to the relation of, 



-i<k,(2) 



-i4>t n (n) > 



,(1,2,.. 



n) 



(38) 



where /^^...(l, 2, . . . , n) is defined by 
A^ 2 ...^„(l,2, ...,n) = 

1(^(1)-^ (2)+ ... - i^MTU^c- (39) 

According to this, AT^ 1 £ 2 ...(1, 2, . . .) is the sum of 
the second-order correlation functions among all time- 
arguments, so that all sequential processes of microscopic 
variables in each term of Eq. H35J) are correlated addi- 
tionally. Actually, the second-order correlation function 
of ((f)£(l)(f)g/ (2))-n RLC is related to the phase fluctuations 
of Eq. (|19|l . One of simple ways for this is to display the 
fluctuation-dissipation theorem in terms of i^ 1 ^ 2 (l,2) 
and then, by comparing it with Eq. 119fl . one obtains; 



+ 

(<^(*)<M*o))hbz, c = 



(Mt)<t>t>(to))o 

0(t-t o )(<M*)<M*o))o 
0(to-*)(&(*)&'(*o)>o 

6(t -t)(Mt)&Ato))o 



(40) 



Sr)(t ) 



S v (t) 



where superscripts of ±=p denote each section of the 
Keldysh space in which both time arguments t and to 
belong. 

The generating functional W is now calculated by sum- 
ming all-connected diagrams in Z. To prevent the diver- 
gence of the average charge at Ao = 0, we perform dia- 
grammatic sum to infinite order, however approximated 
forms in higher-order diagrams cannot be inevitable for 
a simple form of W. Expanding the partition function of 
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Eq. I|35|) and arranging diagrams, we write the approxi- 
mated generating functional as, 



--We^R- 



(41) 



where we omit trivial non-interacting terms. In the first- 
order term, a single diagram contributes the generating 
functional and its self-energy has a form of, 



E«(t,t ) = 5>W°)£{(Mo) 

t 

with a free-environment part of, 

E|(t,t ) = -2ihg d (to,t)b e (t,t ). 



(42) 



(43) 



However, for each higher-order term from the second, 
rich connected diagrams are found. We approximate the 
generating function by taking into account only a single 
diagram in each n_order with its self energies of, 



£(™)(f,t ) = I dld2..dne Kl 

0. ** C 



(i, l)g c (l, 2)£^ (2, 3) flc (3, 4)..E^ (n, t ). (44) 

The idea for such preferred diagrams comes from the 
work of Utsmi et. al.jioj Actually, the partition func- 
tion of Z is the same as that in their work if it were 
not for the electrical environment and time-dependent 
perturbations. So we compose the approximated gener- 
ating functional to recover their results in the absence of 
the electrical environment and time-dependent perturba- 
tions, say, in the case of Ki... — and pt(t) = 0. 

As a result, the generating functional can be written 
in a more compact form as, 



W = -ih \Trllng- 1 } 
= -jftTrtlnG" 1 ] 



71=1 J 



(45) 



where we introduce the full ofield Green's function G 
obeying the Dyson equation; 



G-\t,t Q ) =g-\t,t )-E(t,to) 



(46) 



with a noninteracting Green's function g c (t, to) of the c- 
field. Here, the self-energy E(t, to) is defined by, 



S(t,t ) = £«(t,t ) 
1 fe< 2 ) - E« 9c E« 



(Mo) 



E( 3 ) - hl^g c ^ + isCD^EWffeSW 



E (1) (t,t ) 



(Mo) 



(47) 



which is obtained by comparing Eqs. I|41|l and Q45J1. Since 
Eq. I|4(jll already sum up an infinite series, we calculate 
E(t,t ) to the lowest order of SW(t,t ). This is equiv- 
alent to neglecting additional correlations caused by the 
electrical environment in higher-order terms of the above 
equation. In this case, the self-energy of H(t, to) is ex- 
pressed explicitly with the product of terms representing 
effects of time-dependent perturbations and the electrical 
environment, respectively, as 

H(t,t ) = Er>{t,t )+Zs(t,t ), 
Z e (t,t ) = e^E^toK^H 
£,(t,t ) = e^'^EOfoto) (48) 

and thus £° represents the free self-energy from the en- 
vironment and time-dependent perturbations, 



E°(t,t ) = -2h 2 N ch9d (to,t) | T e | 2 
y^,gki(t,to)g k e(to,t). 



(49) 



IV. 



EXPRESSIONS FOR CURRENTS AND 
NOISES 



For the given generating functional and the self-energy 
of Eqs. (|45J) and l|48(l . we now derive exact expres- 
sions for tunneling currents, averaged charges, and their 
noises. Using the standard procedure, we perform func- 
tional derivatives as specified in Eqs. l2*7|l . I|28[l. J2HJ, 
and H30[l. and transform them into the physical represen- 
tation. The transformations are carried out by adopting 
the Keldysh rotator as, for instance of G, 



G~ 



G 4 
G~ 



1 1 

-1 1 



G A 
G R G K 



1 -1 
1 1 



where superscripts A, R, and K denote its advanced, 
retarded, and Keldysh components, respectively. 

Then, with this rotator the Dyson equation of Eq. I|46|) 
is transformed as; 

G R {t,to) =9?(t,to) 

/OO 
dndng?{t, n)S R (ri, t 2 )G r (t 2 , t Q ) 
-oc 

/oo 
dT 1 dT 2 G R (t,T 1 )E K (T 1 ,T 2 )G A (T 2 ,to) 
-oo 

G A (t 7 t Q ) = G R *(t ,t) (51) 

where g R (t,to) is a retarded component of the free- 
particle Green's function specified in Eq. I|18(l . According 
to these relations, G R ' A has the same causality relation 
as that of g R ' A , and additionally satisfies the sum rule of 
limt-to ^G c (t,t ) = 1 with G c = G R — G A . For G K , 
we show it in terms of G R and G A rather than its integral 
equation by noting that (1 + G R T, R )g^ = G R (g R )~ 1 g^ 
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vanishes. On the other hand, the application of the ro- 
tator to the self-energy of Eq. (|48|) leads to 



are given as, 



Ef (Mo) = Z° e R (t-t )$te K u (*<*°) 

Ef (i, to) = E^ (t-to)^^/^*) (52) 

+ i £° c '(t-to)3e K « + ('■*"), 
E/(t,to) = E^(to,t)*. 

Here, we arrange each term to depend only on the cor- 
relation function of K^ e + (t,t ). Thus, the self-energy of 
E^ is directly related to the phase fluctuation of Eq. l(T$|) 
via Ku + (t,t ) = (MO)iMto -t) - &(0)]) . 

A detailed form of the bare self-energy E° depends on 
a cut-off function for the spectral density of the particle- 
hole propagator be . In this work we use a Lorenzian cut- 
off function p(E) = /(£ 2 + E%) with a bandwidth of 
Eq = Ec as in earlier works. @,[l0] Then, by substituting 
free-particle Green's functions and p(E) into (|49|l . and 
transforming into Fourier forms by 



(53) 



we find, 



EHe) 
E^W 



¥'o(1+-t— ) - Vo(- 



2tt 



2tt 



-i7rcoth 



= E 



OR* 



(e) 



(54) 



where = e— //^ and ^o is a digamma function. Here, Re 
stands for 1/i^ = 4n 2 N ch \T e \ 2 DjD e /R K with energy- 
independent density of states Dj and Di at the metallic 
island and the lead £, respectively. This is tunneling re- 
sistance of the barrier connected to the lead £, so that 
the parallel resistance is given by 1/ Rt = 1 / Rd + V-Rs 
in a = R K /(4ir 2 R T ). 

For static conditions, an analytical solution of the 
Dyson equation in Eq. (|FTJl is easily obtained by trans- 
forming it to the Fourier space. However, for arbitrary 
time-dependent perturbation since the time-translational 
invariance of the self-energy is broken, the solution shows 
up in a series or it is necessary to solve the problem nu- 
merically. One of approximated solutions may be derived 
by noting that £°(i, to) exhibits rapid decaying behavior 
as a function of time interval t — to. If this is the fastest 
time-variation in the problem, the integral equation can 
be approximated to the first-order differential equation. 
Then, the solution becomes similar to that obtained in 
the the wide-band limit. |27j 

With the physical representation of the Green's func- 
tions and the self-energies, results of the functional 
derivatives are summarized as the followings. The av- 
eraged charge in the island and the tunneling currents 



-e~G*(t,i), 



(Q(t))o 

I\{t) = -eUTf(t,t) 
while their fluctuations read 

S Q (t,t ) = ^G K (t,to)G K (t ,t)-(K 



(55) 
(56) 

C), (57) 



S U '(t,to) = e 2 Vi{~«(t,t Q )G K (to,t)6ee> 

+ Kf l ,(t,to)G K {to 1 t)-Tf{t,t )^(to,t)) 
- (K^C) (58) 

where (K — > C) means the change of a Keldysh com- 
ponent to a correlated one, for example, G K — > G c = 
G R — G A . Here, the functions of T and A are defined as, 



rf' c (Mo) 



A«' C (Mo) — 



/ dr 

J — oo 


~?- C (t,T)G A (r,t Q ) 


5?(t,<r)G*' c (r,to) 


3 


p CO 

/ dr 

J — OO 


rf C (t,r)E 


1 a (t, t ) 


If(t,r) 


"l> ( T ^o) 





(59) 



with = — tlf. With these fluctuations, the noise 
spectrum at a frequency u> is defined by, [33, tMJ, |4Cj 

S Q (u) = {{S Q {t,ta)coBu(t-to))) 
Su>(u) = ({Su>(t,to)cosLj(t-t ))) (60) 

where the double bracket denotes the integration of, 



= 31 



lim — 



dir. 



dt 



(61) 



As discussed in the previous section the above results 
obey the charge conservation law even under arbitrary 
time-dependent perturbations, i.e. 



d(Q(t))o 
dt 

d 2 S Q {t,t') 
dtdt' 



= £m 

= ^s, f (M'). 



(62) 
(63) 



This is also a direct consequence of the guage-invariant 
generating functional of Eq. (|45(l . 

The equations of Eq. 1(55 (1 -1(58 ]) are the main results 
of our work together with the generalized self-energies of 
Eq. 148(1 to arbitrary time-dependent perturbations and 
electrical environments. Even though the above equa- 
tions give the exact expressions within the given gener- 
ating functional and provide easier ways in numerical cal- 
culations, their physical meanings are not well revealed. 
So, in the subsequent section we present another forms 
by considering various limits. 
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A. Wide-band limit 

In order to obtain more physically meaningful forms 
we start with defining a spectral function A(,(e,t) as 



At{e,t) 



d T e i < T - t ^ h e- ipe ^G A (T, t) (64) 



and writing Ef^ (e) as, 



Ef (e) = E£(e){2/,( e )-l} 



(65) 



where fe(e) is a Fermi-Dirac distribution function broad- 
ened due to the presence of the dissipative environment, 
otherwise, it is equal to f e (e) = 1/(1 + e^ 6- ^). Then, 
from Eq. I|55(l the average charge can be rewritten as, 



(Q(t))o = e 



J2 fde n e (t,t;e){f+(e)-f e (e)} (66) 



where, by defining a hole-distribution function, //(e) 
1 — /^(e), we emphasize the roles of hole and electron 
contributions. Here, m{t, t;e) is a diagonal value of, 

n e (t,t ;e) = ^e^^^tf (e)A$(e, t)A t (e; t ), (67) 

47T 

which can be interpreted as, owing to its unit, the evolu- 
tion of density of states at the lead I available to occu- 
pying the island. 

As for the tunneling currents, we additionally exploit 
the relation of (1 + E R G R ) = (iWt ~ A )G R and then 
obtain the form of, from Eq. (|56|l . 



— — 5ft / de 

2eR K 



{T[{t,t;e)-T*{t,t;e)} {ft{e)-h{e)} 

-T[(t,t;e){f+(e)-f- e (e)}} (68) 

where the first term in the right-hand side represents 
the particle flux from the lead I and the second is the 
flux from the opposite-side lead £. Here, Tf(t,t;e) and 
T R (t,t;e) are transport coefficients representing tunnel- 
ing of holes (electrons) in the region of positive (negative) 
energies. The function Tf (t, t; e) is a diagonal part of, 

/DO 
dre 1! " ( ' ) E < c (t,r)e-"" (T) 
-oo 

ne{r, t ;e) (69) 

which depends on both sides of the tunneling barriers 
and can be considered as the transmission function from 
the lead I to the other side i. In writing Tf, we adopt 
an approximated form for its further use by substitut- 
ing E^ instead of 2E^. This approximation corresponds 
to the wide-band limit. [27| With this transmission func- 
tion alone, the expression of Eq. I|68|) is consistent with 



the well-known Laudauer formula. |33j However, the ad- 
ditional function Tf gives rise to a deviated form from 
the formula; Tf is given by, 



T f R (t,t ;e) = 8iri 



tft^-Ao ) „,(/./„;,). (70) 



According to this, Tf(t, to; e) is independent of the other- 
side tunneling barrier I contrast to Tf implying not a 
transmission function and, moreover, it does not appear 
in the flux from the other-side lead I. These facts are 
likely to interpret the function as a transport coefficient 
describing the electron or hole flux supplied by the metal- 
lic island. This becomes more apparent if one compares 
T R with Eq. lj6*6"| where the flux by T R is equal to the de- 
crease rate of charges in the island. This flux has no a sta- 
tionary component to be $lT R (t, t; e) = 0, by which our 
expression successfully recovers the Landauer formula in 
a static condition. 

On the other hand, the charge noise is rearranged to 
be, 



Sq(Mo) = '(*>*«) 



Sq (Mo) = 4e 2 / deide 2 n l (t,t \e 1 )n t ,(to,t;e 2 ) 

J • — OC 

{Mei)f+(e 2 )+fn^)fi(^)}, (71) 

which is the sum of all possible electron-hole correlations 
tunneled from both leads. 

The current noise in time-dependent cases is found to 
have a complicated form due to the fluctuations arising 
from various origins and, therefore, we show it under the 
wide-band approximation separating into three contribu- 
tions. The self correlation of the tunneling current at the 
lead I is arranged as, 

5«(t,*o) = S e v (t,t ) + S e r (t,t ) + J^-S%(t,to) (72) 

while the correlation between the different-side tunneling 
currents can be obtained from Eq. I|63l) together with the 
above one. Here, the first term reads, 



S v (t, t ) 



1 



deide 2 L e (t, to] ei)ng{t , t] e 2 ) 



{^(ei)/+(e 2 ) + //(ei)/^2)}(73) 

where Lt(t, to; e) is defined by, 

L t {t, t ; e) = %Tne ipi{t) - ipt{to) e- tt{t ~ ta)/h ^{e). (74) 

Actually, in a static condition this term is proportional to 
the transmission coefficient T F and represents the noise 
arising from backward tunneling after electrons and holes 
starting from different leads, respectively, tunnel into the 
island simultaneously. As a result, the term is not corre- 
lated by the Coulomb interaction in the island, by which 
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it reproduces the results of the co-tunneling theory. | 
The second term of Eq. I|72l) represents the correlation 
between the real tunneling currents; 



si(t, to) 



e 2 R 2 K 



5R 



de\de 2 



T e R (t,t ;e 1 )T[(t ,t;e 2 ){Me 1 )f+(e 2 ) + f+(e 1 )Me 2 )} 
-T R (t,t ;e 1 )T F (t ,t;e 2 ){fe(e 1 )f+(e 2 )+f+(e 1 )f i (e 2 )} 



(Si>,i» - -)T F (t, t ; ei)T/,(t , t; e 2 ) 

{ft(ei)f^(e 2 )+f^(e 1 )fe'(e 2 )}l (75) 



where each transport coefficient of T F and T R represents 
the corresponding current and, thereby, the product of 
them describes their correlation. In a static case, the 
function T R has no a real value, so that this term recov- 
ers the well-known form that represents the correlated 
tunneling processes via the Coulomb interaction. By re- 
calling the term 5y to the self correlation, some of the 
correlations between T R and T F are missing or addition- 
ally incorporated in Eq. (|?5j) . This is a consequence of 
the last term which emphasizes the role of the charge 
fluctuation in the island. In terms of the noise spectrum, 
this term is written as uo 2 Sq(lo). 



B. Adiabatic limit 

For sufficiently small frequencies of driving fields, it 
is expected that the time-dependent fields just vary the 
chemical potentials adiabatically, so that formal expres- 
sions are the same as those in static problems. This is 
really the case if one expands phases in the self-energy as 
Pe(h) = pe(to) +Pe(to){ti -t ) + . .. and neglects higher- 
order terms; 



to; 



n t (t,t ;e) = — Ef(e) \G A (e)\ 2 e^-^, 

T[(t,t ;e) = 4iritf(e)n t (t,to;e), 
L e (t,t ;e) = 87n£f(e)e-^*-*°)/ fi (78) 

and a spectral form of G A {t 1 to) is given by, 



G A {E) 



E-A - X A (E) 



(79) 



Substituting the above expressions into Eqs. (|66fl . 16811 . 
(|71|l . and 172fl . one can find that static results in the work 
of Ref. [T(| are exactly recovered. The function T R has 
no effects on results because its value is imaginary in the 
approximation. 

It is interesting to find the boundary within which the 
adiabatic expressions are valid. To do this, we propose 
the transmission function T F (t,to;e) in Eq. (|69|) as a 
precursor. Then, expanding phases in the self-energy as 
in Eq. i|76|l and requiring the zero_th order term of T F 
much larger than the largest non-adiabatic contribution 
(the first-order term), one obtain a criterion for the valid 
adiabatic expressions as, 



Y,M{t)/R t 



(fru A ) 2 + ±ll 



< i 



(80) 



where the absolute sign means a root-mean-square value 
of a time-dependent function, oja is a frequency of ex- 
ternal perturbation, and a total tunneling rate 70 = 
SE j4 (Ao). According to this criterion, the adiabatic ap- 
proximation is well hold under the case of a smaller ap- 
plied frequency oja than the total tunneling rate 70 to- 
gether with relatively weak amplitudes of perturbations, 
as usually expected. However, in general, it also depends 
on temperature, capacitances of the tunneling barriers, 
as well as the excitation energy A . 



miM) = e imtl) !Mti,i 2 )e-^) 



(76) 



Then, this changes simply the chemical potential in Eq. 

as; 

-> to(to) = $-ta$(to)-(5t, D -C D /C E )eVg (77) 

where we separate the potential difference across a tun- 
neling barrier I into an alternating part of v\(t) and a 
direct part with a static external voltage V§ of Vu(t). 
Here, we use the time argument to to clarify the fact that 
the chemical potentials in noise calculations depend on 
the reference time to rather than another time argument 
t. 

Under the adiabatic approximation, the effective den- 
sity of states and the transport coefficients are simplified 



C. Orthodox results 

Under the adiabatic regime we further approximate the 
formalism by assuming very opaque tunneling barriers. 
Then, | G A (e) | 2 behaves a nearly (5-function, reflecting a 
long life time of charged states. This makes possible the 
integration in the formalism as, 

/oo 
dcF(e)\G A (e)\ 2 =-F(A ) 
-00 7o 

/OO 
deF(e)\G A (e)f\G A (e-u,)\^ 
-OO 

■[F(A ) + F(A +^)] (81) 



7o(4 7 2 +^2) 



where we neglect a renormalization effect on the excita- 
tion energy A . As a consequence, results become iden- 
tical to those based on the orthodox theory. 0, ^(j 
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V. APPLICATIONS OF FORMALISM 

As applications of our formalism, we now examine the 
performance of rf SETs numerically. For this, we consider 
detailed geometry of the circuit which is characterized by 
the impedances [Zq and Zp>) and the voltage sources (Vb 
and Vq) in Fig. According to measuring processes of 
signals, there are two kinds of SETs; reflected and trans- 
mitted types. In this work we assume the reflected type 
of rf SETs where a detector measures reflected signals 
from a SET. 0, 0, 0] Applications to the transmitted 
type0, ^3 is basically identical with slightly modified 
electrical environments. 

Then, an equivalent circuit for the rf SET is given by 
writing the impedances as 



Zq(w) = 

Z D (ui) = Rq- - 

\ — uj z + lui/Qi 



(82) 



where i? is the coaxial cable impedance (typically 50O), 
w = uj/lur is a normalized frequency with a resonant 
frequency ujr of the tank circuit, and Qp is its quality 
factor. We model microwaves propagating the coaxial 
cable as a sinusoidal form, v~~ (t) — Vi n cos(uJAt) with am- 
plitude and angular frequency Vi n and uj a , respectively, 
and also consider a static voltage provided by a bias 
tee. Then, the equivalent voltage source of Vd (t) is given 
by, in its Fourier component, 



VbM = 



Zd(uj) 2v (u>) 



6(lv)V° 



(83) 



R 1 + iujQf 
or in a real-time space, 

V D {t) = 2r 1 o(uj A )v in cos(LU A t + cp(uj A )) + Vg (84) 



with the phase and amplitudes defined by 

1 



rj (u)e 



1 — uj 2 + iCo/Qp 



(85) 



A gate voltage Va(t) represents a signal to be measured, 
for instance, the voltage induced by time-varied charges 
of qubits. If the signal is slowly varied in time, it can be 
treated in the adiabatic way like the change of a static 
voltage Vq , which in turn gives rise to the modulation of 
the excitation energy Ao adiabatically. In the followings 
we assume this case where Vc(t) is no longer a voltage 
source, i.e.,Va{t) = 0, but a parameter for the excitation 
energy Ao. However, in general it should be treated time- 
dependent fields for considerations of gate charges varied 
rapidly. 

From the voltage sources, the potential differences 
across the tunneling barriers from Eq. (|13|l become 



v b D (t) 



1 




(86) 



with Va{t) = 0. Here, the induced voltages from the 
tunneling currents are calculated from Eq. I|llfl as, 



vh(t) 



1 

Ta 



Ta 



drlUr) [R + 2^Z D {ij A )e^^ 



1 f A e l 
V G (t) = — / dr[J|,(r)+^(r)]29f 7 



(87) 



iuj A C G 

where, due to resonant properties of the tank circuit, 
the lowest two harmonics of the tunneling currents are 
taken into account and T A = 2tt/u) A is a period of the 
microwave. 

As for output signals, we consider a reflected voltage 
v+(t) from the SET; 

v + (t) = v in cos[u A t - 2(p(uj A )} - R I 
+ X cos(uj A t) + Y sm(uj A t) 
+ higher harmonics. (88) 

Here, we separate a pure reflected component of the mi- 
crowave (the first term in the right-hand side) from its 
responses to the single-electron transistor. Each term de- 
pending on Iq, X , and Y is originated from the tunneling 
current, 



I 1 it) 



1 



Cv) lD{t) c s 

Io+Ii cos(u) A t) + I[ s'm(LJ A t) 
higher harmonics, 



(89) 



and their Fourier components are related to each other 



as, 



X 
Y 



Rq 



VP 



(90) 



To model a detector we consider the average of observ- 
ables multiplied by cos tot or sin Lot over time. 0> 03 Es- 
pecially, we focus on a homodyne detector measuring the 
amplitude X obtained from the reflected voltage v + (t). 
As indicated in Ref. [42|. since the amplitude X is usu- 
ally much larger than Y for Qp 3> 1 implying a small 
reactance of tunneling barriers, it is a good approxima- 
tion to express the performance of the rf SET only in 
terms of X. After many numerical simulations we find 
that this is also the case for our system. 

Then, the noise associated with X is derived as, 



2R 2 

Sxx(u) = j- ~ 2 \2 i ~2 ir>2 S i( u >) 

(1 - uj z ) z + bJ z /Q F 



(91) 



with noise of a tunneling current Si{uj) 

Sj(uj) = £(4 

w 

((Su'(t + to, to)[cos{2uito + 2lp(lj)} + cosoj 



Cd v x C d 
\D--pr){n',D--^r) 
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It is noted that the current noise of Si (to) is deter- 
mined from the average weighted with a factor cos{2ujt + 
2ip(u)} as well as cos wt. This is a consequence of model- 
ing the homodyne detector. In static cases, the average 
with the former gives zero because the current fluctua- 
tion of Sw(t + to, to) is invariant under time-translation, 
and then Si just represents noise of the total current 7* 
of Eq. (|89|) . However, in general since the invariance is 
no longer valid under time-dependent conditions as easily 
checked in Eq. I|72[l. the average has a finite value and, as 
a result, the noise of X contains additional contributions 
at every multiples of the half- frequency to a/ 2. 

The zero-frequency approximations correspond to set- 
ting of to = in cos cot for considering a small frequency 
lu and then our expression of Eq. I|92[l becomes similar 
to Eq. (30) of Ref. 42]. 

Our numerical examinations are fulfilled under several 
limitations due to the two-states approximation. So, the 
range of parameters for operating points such as the ex- 
citation energies A (or gate charges qo), the amplitudes 
Vi n and frequencies to a of microwaves, and the DC bias 
voltage Vq should be restricted not to occupy higher or 
lower charged states. For this we consider the excitation 
energies in the range of — Ec < A < Ec (or < qo < e) 
and microwave energy Hu>a much less than Ec- When 
the frequency of a microwave is tunned to be resonant to 
the tank circuit (i.e., to a = ^r, which is also assumed for 
our numerical calculations), the maximum amplitude of 
2QFVin is delivered to the SET. Then, not to excite other 
charged states, the applied amplitudes should satisfy the 
inequality of 



eV D \+2QpeVin < min 



cv 



2E c (l 



, D <-y£ — 

<Zo_ 1 
e 2 



(93) 



which can be derived from simple electrostatic consider- 
ation. According to this, energy provided by the drain 
voltages (the left hand term in the above equation) is 
restricted to be approximately less than Ec and 2Ec at 
completely blockade points (qo = 0, e) and a degenerate 
point (qo = e/2), respectively. 

Hereafter we use local units for calculated results and 
system parameters. We display all quantities in units 
of charging energy Ec, but in units of an ohm for re- 
sistances. In other words, thermal energy fcgT, excita- 
tion energy Ao, amplitude evi n and frequency Hloa of 
microwaves, and DC bias voltages eVj^ are measured in 
units of Ec as well as currents are measured in units 
of eEc/h and capacitances in units of e 2 /Ec- In these 
units, the total capacitance becomes Cs = 0.5. So, 
the situation such as Cd = Cs = 0-2 in the fallowings 
means a small gate capacitance and, therefore, gives rise 
to somewhat symmetric distribution of potential differ- 
ences across the tunneling barrier from Eq. I|8t)|) while 
the case of Cd = Cs — 0.01 represents asymmetric 
distribution. Our numerical calculations are performed 



for i?o = 50f2, a zero static voltage Vjj = 0, and a 
temperature of fc^T = 0.027?c (unless mentioned oth- 
erwise) . Even though the dependence of the SET perfor- 
mance on t emp erature and finite static voltages are also 
interesting, |l5|,|42j we omit it for simplicity. 



A. Environmental effects 

Firstly we discuss effects of the electrical environments. 
For the given electric environment of Eq. <|82|l . its effects 
are manifested into the self-energy of T,g(E). In general, 
the presence of electrical environments broaden the self- 
energies and one can identify this via the Fermi-Dirac 
distribution functions of Eq. I|65l) . In Fig. Q , we show 
the Fermi-Dirac distribution functions for two cases of 
resonance frequencies for various quality factors, in which 
the resonance frequency is comparable to thermal energy 
in (a) and much larger than it in (b). Characteristic 
behavior is that the distribution function becomes more 
depleted around the chemical potential (/u° = 0) as the 
quality factors and the resonance frequencies increase. 
This is resulted from the energy-emitting and absorbing 
spectrum of the environment. As indicated in Ref. |25| a 
small quality factor means rather rapid damping in elec- 
trical dynamics of the environment or Ohmic behavior, in 
which the spectrum has a peak around a zero irrelevant of 
the resonance frequency. Whereas, a large quality factor 
leads to the environment with a single mode case where 
quantized energy equal to hujR is incorporated into the 
spectrum and produce peaks at every multiple of TujJr. 
Thus, in the case of the large quality factor, tunneling 
is mediated to the environment by emitting or absorbing 
energy quanta of huiR. As a result, if thermal energy ksT 
is less than it and cannot excite the environment, the en- 
ergy quanta should be provided externally, which in turn 
reduces the effective number of particles for tunneling. 

To emphasize the effect of the depletion in the particle 
distributions we show calculated tunneling currents for 
static cases in Fig. 0}(c) using parameters in (6). For a 
small quality factor (Qf — 10), tunneling currents are 
very similar to that for the free environment. However, 
for a large value of Qf = 80, noticeable differences are 
found. This difference becomes more enhanced for larger 
resonant frequencies over the thermal energy and more 
asymmetric geometry of tunneling barriers like small ca- 
pacitances of Cd and C's compared to Cq. 



B. Effects of photon-assisted tunneling 

Next we discuss effects of time-dependent perturba- 
tions on transport by comparing results calculated from 
the exact and adiabatic expressions. For this we con- 
sider sufficiently large amplitudes and high frequencies 
of a microwave compared to a temperature as well as 
asymmetric geometry of barriers. Calculations are per- 
formed for two cases of tunneling resistances; one is large 
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FIG. 4: We show the Fermi-Dirac distributions at the drain 
for two electrical environments. In (a), we choose a resonant 
frequency of Hujr — 0.04 at a temperature of fcsT = 0.02 while 
in (b) hujR = 0.08 and fcsT = 0.02. For various quality fac- 
tors of Qf = 10 (dotted line), 50 (dashed), 80 (dot-dashed), 
results are compared with that in the case of the free environ- 
ment (solid). In (c) DC tunneling currents are calculated as 
a function of DC drain voltage Vd with parameters in (b) at 
qo = e/2. Here we use Rd = Rs = 13.1fcS7 (ao = 0.1) under 
a nearly symmetric configuration of barriers, Cd = Cs = 0.2. 

resistance of Rd — Rs — 655fcSl (ao = 0.002) and the 
other is a small one, 13. (ao = 0.1). So, the case of 
the former is believed to be dominated by sequential tun- 
neling in transport while with the small resistance both 
co-tunneling and sequential processes are expected to be 
contributed. 

Fig. [S] shows calculated results of the currents and 
noises for the above two cases as a function of gate charge; 
for the smaller (large) in the left (right) column. In each 



a =0.002 a„=0.1 




FIG. 5: We plot the first harmonics of a tunneling current 
(in units of Ec/eRr), current and charge noises (in units of 
Ec/Rt and e 4 Rt / Ec , respectively), and the criterion func- 
tion as a function of gate charge qo for two sets of param- 
eters. In the left column, tunneling resistances are chosen 
to be Rd = Rs = 655fcf2 (ao = 0.002) and a frequency of 
a microwave Tiuja = fiujR — 0.01 while in the right column 
Rd = Rs = 13.1fen (ao = 0.1) and Hcja = 0.2 are used. 
We compare results calculated based on the exact, adiabatic, 
and orthodox formula by displaying them in solid, dotted and 
dashed lines, respectively. Other parameters are hsT = 0.01, 
C D = Cs = 0.1, Q F = 6, and 2Q F ev in = 0.5. 



column we additionally distinguish results depending on 
the formalism used; the exact (solid), adiabatic (dotted), 
orthodox (dashed) formalism, respectively. Thus, the 
differences between the exact and adiabatic results rep- 
resent the existence of photon-assisted tunneling while 
those between adiabatic and orthodox results just em- 
phasize effects of co-tunneling in the classical limit (for 
the small resistance the orthodox theory is known to be 
invalid, nevertheless we also show corresponding results 
just for comparison). 

From the calculated results, one can see that the sys- 
tem is easily driven into the non-adiabatic regime around 
the degeneracy point go = e/2 in both cases of tunnel- 
ing resistances implying the fact that tunneling occurs 
in photon-assisted ways and, there, tunneling currents as 
well as noises are largely enhanced relative to the adi- 
abatic results. This tendency is more apparent in the 
case of the large resistance. However, apart far from 
the degeneracy point, effects of time-dependent pertur- 
bations are reduced to classical ones to give nearly identi- 
cal results independently of the formalism used. In other 
words, the strength of the photon-assisted tunneling be- 



15 



1.0 
0.8 

0.6 
0.4 
0.2 
0.0. 



- 


i i mi 






















(a) 


///// 


— ^-v/ 





a„=0.002 



a„=0.1 



-0.5 



0.0 



0.5 

E 



1.0 



1.5 




FIG. 6: We plot the transmission function T^(t, t; E) at the 
drain with parameters used in Fig. ©-(b). For gate charges 
of qo — 0.2e and 0.5e, its contour plot is depicted in (a) and 
(b), respectively. Here, a time t is measured in units of a 
period of a microwave Ta- 
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FIG. 7: Under the same condition as those of Fig. |S] we 
plot the frequency dependence of current (top) and charge 
(bottom) noises for ao = 0.002 (left) and 0.1 (right panel), 
respectively. In each figure, results by the exact, adiabatic, 
and orthodox formula are represented by solid, dotted, and 
dashed lines, respectively, at a gate charge qo — 0.4 (cur- 
rent and charge noises are measured in units of Ec/Rt and 
e 4 Rr I Ec , respectively). Sharp peaks correspond to an infi- 
nite averaging time T in Eq. (1611 . otherwise smeared ones 
may be obtained. 



comes weak far from the degeneracy point. 

One of possible explanations for the different strength 
of the photon-assisted tunneling depending on gate 
charges qo and tunneling resistances is the resolution of 
photon energies Hloa seen by particles. As inferred from 
Eq. I|79l) since particles in the island is decayed with a 
rate of 70 = 3£ A (A ), its dwelling time can be regarded 
to be inversely proportional to the rate. Along this as- 
pect, the left-hand term in Eq. I|8l)|l can be interpreted as 
rough estimator of the resolution for energy quanta YvjJa- 
We plot it in the bottom panel of Fig. [SJfor both cases of 
the tunneling resistances. According to the figure, when 
the value is less than about 0.1, the calculations show 
nearly identical behavior between the adiabatic and ex- 
act results implying poor resolution for energy quanta. 
This fact is also confirmed by studying the transmission 
function T F of Eq. lff]9")l . Fig. El is a contour plot of 
the transmission function in the energy-time space for 
two different points of gate charges at which the reso- 
lution is predicted to be poor (at go = 0.2e) and good 
(at qo = 0.5e) in Fig. El(b). As shown in the figure, 
their time and energy dependences are strikingly com- 
plex. However, one can recognize the distinct feature 
different from each other. In (a) for a given time the 
height of the transmission function is changed continu- 
ously as a function of energy while in (b) it exhibits rather 
discrete behavior. Actually in the case of (b) the distance 
between successive maximums is equal to photon energy 
of foui a, which is well resolved enough to see photon-side 
bands. 

For the frequency dependences of the current and 



charge noises, we examine them at a gate charge qo = 
0.4e of Fig. Eland show results in Fig. 0in the range 
of low frequencies. In the current noises, the differences 
at a frequency uj = to a between the adiabatic and exact 
results are found to be retained in the whole range of 
frequencies accompanying some structured behavior in 
both tunneling resistances. A straightforward explana- 
tion for these differences is complicated by various corre- 
lations among transport coefficients of ni, Tf , T^, and 
Lg in Eq. l|72Jl where they all are found to play roles 
somewhat importantly. As for peaks located at multi- 
ple of a frequency ui a I '2, they are apparently arisen from 
time-dependent properties of external perturbations as 
noted in Eq. (|92|l . Usually, the largest peak is found at 
w = uja, however, in the exact results peaks at u> = ui a 1 2 
and 3lua/2 are also appreciable contrary to the adiabatic 
cases. This means that in photon-assisted ways the cor- 
relations between different harmonics of the above coef- 
ficients remain large while they are unimportant in the 
adiabatic limit. In the bottom of the figure, we show 
the charge noises as a function of frequency. Due to 
a large frequency scale we do not distinguish their dif- 
ferences clearly. However, in term of lj 2 Sq(ui), we find 
that each result shows an appreciable deviation from the 
others. Actually, the function of lo 2 Sq(oj) is related to 
the current noises because the last term of Eq. I|72|l is 
approximately equal to it neglecting geometrical factors 
between capacitances. By comparing it with the current 
noises, we find that a large portion of the difference in the 
current noise can be associated with those in the charge 
noises. 
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For high frequencies, the noises are found to be nearly 
identical to equilibrium noises calculated with no exter- 
nal perturbations. Namely, the Johnson-Nyquist noises 
are dominant. 



C. Charge sensitivity 

The sensitivity to gate charges is one of estimators for 
the performance of a rf SET as an electrometer. If one 
measures the amplitude X assuming a homo dyne detec- 
tor, the charge sensitivity is calculated as, @,E2 



Sq = 



dXj dqo 



dh/dqo 



(94) 



showing the dependence on the current noises Si and 
the response function dli/dqo. In Fig. 0we plot the 
sensitivity in the qo-Vi n plane using the three different 
formula. Aiming at the simulation of Ref. |l5j . we use 
similar system parameters to the experiment except for 
Cd = Cs = 0.1 (unfortunately not specified in the ex- 
periment). To clarify calculated results we omit the re- 
gions of poor sensitivity, which correspond to smaller rf- 
wave amplitudes than the Coulomb blockade thresholds 
(lower left corner) and much higher rf-wave amplitude 
over the threshold with small excitation energies (upper 
right corner). By comparing Fig. |U-(a) and (b), one can 
see effects of photon-assistant tunneling to the sensitiv- 
ity. That is, the region of good sensitivity (for example, 
within Sq = 2.5) is predicted to slowly vary as a func- 
tion of rf-wave amplitude than in the adiabatic limits, 
while it exhibits a rather narrower region as a function 
of gate charge. Calculated optimum sensitivities (mini- 
mum value of Sq) are also found to have different values, 
Sq = 2.12,2.38, and 1.83 for the exact, adiabatic, and 
orthodox formalism, respectively, with a nearly same op- 
erating point of (q ,2Q F Vin) = (0.47,0.12 - 0.17). The 
slightly better optimum-sensitivity in Fig. [SJ( a ) com- 
pared with that in (b) comes from the enhanced value of 
the response dI q /dqo by photon-assisted tunneling. How- 
ever, photon-assisted tunneling does not always give the 
better sensitivity because it also enhances current noises. 
We find that its role for the sensitivity depends on sys- 
tem parameters. In the orthodox result, the sensitivity is 
always predicted to be better than results of the others 
because the absence of co-tunneling gives larger values of 
the response dI q /dqo. 

Using the experimental parameter of Eq = 178/j.eV, 
the calculated optimum sensitivity Sq — 2.12 corresponds 
to 4.1/ie/VHz, which is lower than the measured value 
of 47/ie/VHz by an order. This discrepancy between the 
theory and experiment may be attributed to the pream- 
plifier noise and local heating of a SET.|42l| 

In Fig. |21 we show the dependence of the optimum 
sensitivity on tunneling resistances for different applied 
frequencies of microwaves and quality factors. Firstly, 
Figs. E}(a) and (b) are results for Qp — 6 and 50, re- 
spectively. By comparing different kinds of lines in each 
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FIG. 8: Contour plots of the charge sensitivity Sq (in units 
of e^/h/Ec) in the space of rf amplitudes and gate charges 
are shown. The used parameters are huj = 0.04Ec, Qf = 6, 
and Rd = Rs = WOkQ at a temperature of ksT — 0.02Ec- 



figure which means different frequencies of microwaves at 
a resonant condition, one can see that larger frequencies 
give rise to worse sensitives in a wide range of tunneling 
resistance. Since a large driving frequency corresponds 
to better-resolved energies of external fields, associated 
photon-assisted tunneling is expected to decrease the sen- 
sitivity due to the enhanced current noises. On the other 
hand, for a given frequency, the sensitivity as a function 
of tunneling resistance is found to show different behavior 
depending on the quality factor; monotonically improved 
results are obtained for Qp = 6 as tunneling resistance 
decreases while there is optimal resistance for the best 
sensitivity for Qp = 50. For easier understanding of our 
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FIG. 9: The optimum charge sensitivities are plotted as 
a function of tunneling resistance for two quality factors of 
Qf = 6 in (a) and 50 in (&), respectively, together with dif- 
ferent frequencies of a microwave, ftuu =0.1 (dotted), 0.04 
(solid), and 0.005 (dashed). The dot-dashed lines denote the 
analytical result of Eq. 1951 . Here, we assume Rs = Rd, 
Cd = C s = 0.1, and k B T = 0.02. In (c) we re-plot the 
result of (b) represented by the solid line, and examine its 
change without each calculational factor one by one; adia- 
batic approximation (dotted line), no electrical environment 
(square), symmetric geometry of Cd = Cs = 0.24 (dashed 
line), and calculation neglecting a self-consistency (triangle). 

results, an analytic expression from Ref. |42| . 

Sq ~ 2Me^J{R D +R s )C^ -|- (95) 

is also plotted with dot-dashed lines in the figure even 
though it gives difference values from ours due to the or- 
thodox theory. According to this expression, since the 



— 1/2 

sensitivity is proportional to a ' for Rp = Rs, one 
can see that calculated results are approximately scaled 
as the same power law. However, for Qp = 50, the power 
law is no longer hold and for small tunneling resistances it 

1/2 

is scaled as even a . This behavior can be understood 
through a simple circuital analysis as in Ref. ^3] By 
replacing the tunneling barriers by resistors, it is found 
that the best sensitivity is achieved for series resistance 
equal to Q f Rq at which input impedance of microwaves 
is matched to that of the tunneling barriers. Thus, for 
Qf = 6 the matching condition occurs at ao = 1.5 while 
for Qp — 50 it is expected to be ao = 0.02. As shown in 
Fig. this condition well agrees with the calculated 

result for the small frequency even though for high fre- 
quencies it is slightly deviated due to the non-adiabatic 
effects of microwaves. 

In Fig. EJ- (c) we show the dependence of the charge 
sensitivity on tunneling resistances by omitting a calcu- 
lational factor one by one to emphasize its role. Accord- 
ing to results, considerations of self-consistent and non- 
adiabatic schemes are found to be crucial for sensitivity 
calculations while effects of the environment is relatively 
unimportant. In addition, we find that symmetric geom- 
etry (dashed line) benefits the charge sensitivity in the 
whole range of tunneling resistance. 

VI. SUMMARY 

In this work, we develop a formalism for a radio- 
frequency single-electron transistor taking into account 
electrical environment, higher-order co-tunneling, and ar- 
bitrary time-dependent perturbations. Assuming large 
charging energy, we use a two-charged-state model in 
a metallic island and solved the problem based on the 
Schwinger-Keldysh approach combined with a generat- 
ing functional method. We calculate an approximated 
generating functional by summing diagrams in infinite 
order and give exact expressions for current, charges in 
the island, and their noises within the generating func- 
tional. By defining generalized transport coefficients, we 
write the derived expressions in terms of them, and show 
that tunneling currents in time-dependence cases have a 
generalized form of the well-known Landauer formula. 

As application of our formalism, we examine tunneling 
currents, its noises, and the charge sensitivity of rf-SETs 
by accounting for a detailed tank circuit. Firstly, effects 
of the electrical environments are found to be relatively 
small, as expected, in cases of microwaves delivered via 
an coaxial cable with impedance 50f2. However, for a 
large quality factor and a large resonant frequency its ef- 
fects become large and cannot be ignored. Secondly, ef- 
fects of photon-assisted tunneling are manifested to both 
enhanced responses and noises of rf SETs. However, due 
to the larger enhancement of the noises, photon-assisted 
tunneling is not helpful to the charge sensitivity. As a 
consequence, with experimental parameters of Ref. |15| , 
we obtain the charge sensitivity of Sq = 4.1^e/\/Hz, 
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which is larger than that in the orthodox result, however, 
still much smaller than the measure value of 47//e/VHz. 
Finally, we discuss the charge sensitivity depending on 
various sets of parameters. Especially, we focus on its 
change as a function of tunneling resistance, and find that 
the charge sensitivity for small quality factors is scaled 



like 



,-V2 



as in the analytic result proposed by the pre- 



vious work. Whereas, for large quality factors the power 

1/2 

law is no longer valid and it is proportional to even a Q 
in the range of small tunneling resistance to show optimal 
resistance for the best sensitivity. 
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